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ABSTRACT 

We test anisotropic dark energy models with the 7-year WMAP temperature observations data. 
In the presence of imperfect sources, due to large-scale gradients or anisotropics in the dark energy 
field, the CMB sky will be distorted anisotropically on its way to us by the ISW effect. The signal 
covariance matrix then becomes nondiagonal for small multipoles, but at t > 20 the anisotropy is 
negligible. We parametrize possible violations of rotational invariancc in the late universe by the 
magnitude of a post-Fricdmannian deviation from isotropy and its scale dependence. This allows to 
obtain hints on possible imperfect nature of dark energy and the large-angle anomalous features in the 
CMB. A robust statistical analysis, subjected to various tests and consistency checks, is performed 
to compare the predicted correlations with those obtained from the satellite-measured CMB full sky 
maps. The preferred axis point towards (7,6) = (168°, —31°) and the amplitude of the anisotropy is 
voq = (0.51 ± 0.94) (la deviation quoted). The best-fit model has a steep blue anisotropic spectrum 
(n de = 3.1±1.5). 

Subject headings: cosmic microwave background — cosmology: observations — methods: numerical 



1. INTRODUCTION 

In the past two decades great advances have been 
made in observational cosmology. The most striking sin- 
gle discovery is the present acceleration of the universe 
expansion, now confirmed by many independent exper- 
iments. The most powerful probe of precision cosmol- 
ogy is t he observations of the cosmic microwave bac k- 
ground (jBennett et all (j2003l ): iHinshaw etall (|2007f) ). 
which seem to support the model of universe which at 
large scales is flat, isotropic and homogeneous, as firmly 
predicted by inflation. However, at more subtle level 
there seems to be also hints of substantial anisotropy. 
Such would imply violation of the cosmological principle, 
perhaps as striking change of paradigm as the introduc- 
tion of dark energy. Presently the evidence for anisotropy 
is debatable, but the bounds can definitely be expected to 
improve with the Planck experiment. Therefore it is ex- 
tremely interesting to study theoretical links between the 
acceleration and anisotropics, in part icular, the possibi l- 
ity to constrain them observationally iCopi et aLl (|2010f) . 

Several distinct statically anisotropic features have 
been reported in the data analysis of the CMB sky. 
Among the most curious is the p resence of hemispherical 
asymmetry (jEriksen et al.l ((2004)). Recent investigations 
exploiting the five-year WMAP data have found that the 
evidence for this asymmetry is increasing and extends to 
much smaller a ngula r scales than previous ly believed to 
([Hansen et all (pOOct ): iHoftuft et all ([20091 )1 Alignment 
of th e quadrupole and qctup ole, the so called Axis of 
Evil (jLand fc Magueiiol (j2005D ) could also seem an un- 
likely result of statistically isotropic perturbations, even 
without taking into account that these multipoles happen 
also to be aligned to some extent with the dipole and with 
the equinox. In the CMB spectrum, the angular correla- 
tion spectrum seems to be lacking power at the largest 
scales. The alignments seem to be statistically indepen- 



dentof the the lack of angular po wer (|Rakic fc Schward 
(120071)). For other studies, see (jPrunet et all (|2005D : 



IGordon et all ((20051) ). 

It is natural to associate the apparent statistical 
anisotropy with dark energy, since the anomalies occur 
at the largest scales, and these enter inside the horizon 
at the same epoch that the dark energy dominance be- 
gins. The paramount characteristic of dark energy is its 
negative pressure. One may then contemplate whether 
this pressure might vary with the direction. Then also 
the universal acceleration becomes anisotropic, and one 
would indeed see otherwise unexpected effects. These 
would presumably be strongest at the smallest multipoles 
of the CMB since they describe the large angular scales 
which are most directly affected during the late epochs 
of the universe. Specifically, as the photons travel from 
the last scattering surface towards us, the their temper- 
ature gets blue- and rcdshifted as they fall in and climb 
out of the gravitational wells, respectively. When the 
potentials evolve, there is a net effect in the temperature 
of the photons: this is the integrated Sachs-Wolfe effect 
(ISW). Furthermore, if the average evolution of the po- 
tentials was not the same in different directions of the 
sky, the effect would be anisotropic. However, to explain 
the lack of large-angle correlations, there should occur a 
cancellation with the Sachs- Wolfe effect from the poten- 
tials last scattering surface that typically contribute to 
the l arge angles with simi lar order of magnitude as the 
ISW IAfshordi et al.l (pOOflh . 

The potentials parameterising the perturbations of the 
metric, can be written in the longitudinal gauge as 



ds 2 



a 2 (v) [-(1 + H ] )dif + (1 - 2(j))dx i dx l 



(1) 



The Poisson equation relates the spacetime curvature <j) 
to the matter sources. As is well known, in the ab- 
sence of anisotropic stress the potentials 4> an d ip are 
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equal. Thus, detection of inequality of these potentials 
in the present universe would indicate the presence of 
imperfect energy source, eithe r in form of dark energ y 
fluid or modificatio n of gravity (IKoivisto fc Motal (120061): 
iDaniel et all (|2008| ) : IManera fc Motal (|2006| ): lDvali et al.1 
(|2000h ). Clearly, the difference of the potentials can be 
constrained muc h tighter at Solar system than at cos - 
mological scales (iMota et all (1200 7]): IDaniel et al.1 (I2009D; 
iFerreha fe Skordisl (I2010f ) : ISkordisI (l2009D : lHu fe Sa"wickl 
((20071 ) ) . In the present study, we consider the possibil- 
ity that the relation of the anisotropy described by the 
difference of the potentials does not cancel out on the 
average, i.e. that the anisotropy is statistical. 

This amounts to promoting each Fourier mode of 
the potentials to depend not only on the length but 
also on the direction of the wavevector. This is a 
generic prediction for perturbations in a non-FRW uni- 
verse and al so for non-scalar field models, in particular 
vector fi elds lArmendariz-Piconl (I2004D ; IKoivisto fe Motal 
(|2008aD ; IZu"ntz et al.l (|2010l ); lLi et al.f(|2008f) . It has been 
considered if the cosmological fine-tunings could be more 
naturally alleviated with a dynamical dark energy com- 
ponent when this is modelled with a more general field 
than a scalar. Vect or fields dynam i cs could accelerate 
the un iverse today (|Kiselevl (|2004D ; Uimenez fc Marotol 
20081) ) ha v ing ph a ntom evolution withou t UV pathology 



Rubakovl (|2006fk iLibanov et all (120071 ) ) possibly con- 



necting the acceleration wi th the electromagnetic scale 
(jjimenez fc Marotol (|2009al lbT) : Uimenez et all (|2009h V 

There has also been i nterest o n anisotropies in i nflatio n 
(IGumrukcuoglu et all |2006l) : lAckerman et al.l (|2007D ; 
IBoehmer fc Motal (120081)17 and their c omparison with the 
data ( Groeneboom fc Erikse n (2008) ;lGroenebqom et all 



(|20l0h : lArmendariz-Picon fe Pekowskvl (|2008l )). The 
anisotropy in the primordi al spectrum co u ld b e 
generated by v ec tor fie lds (jGolovnev et all (|2008D ; 
IKoivisto fc Motal (l2008b|)) o r more general n-form s 
(iGermani fc Kehagiasl (120091) ; IKoivisto et all (|2009l ); 
IKoivisto fc Nunesl (|200S )). In particular this can be 
robustly realized b y the vector c urvaton paradigm 
Dimopoulos et all (|2010l )) which IDimopoulos et all 



20111 ) recently implemented within D-brane inflation 



in type II string theory by taking into account the 
U(l) gauge field that lives on the brane. Perturba- 
tions have be en also studied in a n isotropically inflating 
backg rounds (jPereira et all (|2007D ; IGumrukcuoglu et all 
(120071)) and in the s h ear-free c osmo logies considered in 
IKoivisto et"aLl (|2011[ ); IZlosnikl (|2011l ). where the expan- 
sion is isotropic but the spatial curvature depends on the 
direction. These homogeneous but anisotropic universes 
coul d emerge by tunneling f rom a lower-dimensio nal vac- 
uum lAdamekeian (|2010f) ; IGraham et all (|2010D . 

In the presence of such variety of possibilities, we 
choose to rather employ a general parametrisation than 
study a particular model. To that purpose, we parame- 
terize directly the angular variation of the gravitational 
potentials. This can be seen as a step towards a more 
complete anisotropic post-Fricdmannian parametrisation 
of the deviations from standard GR ACDM cosmology, 
inspired by the recent development of a fully consis- 
tent par ametrisation encompassin g statistically i sotrop ic 
models (jFerreira fc Skordisl (|20Tot) ; iBaker et all (|201ll) ). 
In Section [2] we derive the signal covariance matrix in 
the presence of generalised perturbation sources. In Sec- 



tion [3] we describe our parametrisation of such sources 
and their interpretation as a as anisotropics of the dark 
energy field or as some spontaneous anisotropisation of 
the CMB radiation. In Section|3]we discuss the analysis, 
and the method we use is described in detail in Section 
[H and finally, the results are presented in Section [5] 

2. CMB FROM ANISOTROPIC SCALAR SOURCES 

The temperature anisotropy field is conventionally ex- 
panded in terms of the spherical harmonics and on the 
other hand considered in the Fourier space 

6( X ,e,r ? )=^ E a *» Y tr»= / 7^e jkx ,5(k)e(k,e,7 ? ), 

(2) 

where we have normalized the transfer function 0(k, e, r/) 
with respect to the initial amplitude <5(k). One 
makes contact between the two expansions by using the 
Rayleigh formula 

oo 

e ,x-k = J2 % \2t + l)j £ (kx)Y; m (k)Y em (x), (3) 

£=0 

together with the addition theorem for the spherical har- 
monics 

4-7T ^ 

Pi{k ■ P) = ^— £ Y; m (k)Y t , m ,{p) (4) 

m=-£ 

for the second equality in Eq.Q and then, by exploiting 
the orthonormality of the spherical harmonics, picks up 
the coefficients in the first equality in Eq.@. These are 



a hn — l 



d 3 k 
2^2' 



<5(k)r/ m (fc)e ( (k). 



where we have defined 



e,(k)= / j/(*rfo))e(M)d»7. 



(5) 



(6) 



We assume, as usual, that the primordial spectrum of 
perturbations is statistically isotropic, 



(S(k)S*{k')) =P(fc)(27r) 3 <5 3 (k-k'). 



(7) 



However, we allow the transfer function an anisotropic 
part, 

Ql(k) = e i (k)+w(k-h)ef-(k). (8) 

The second term can then incorporate the anisotropic 
ISW contribution from dark energy. We are then inter- 
ested in the correlators 



2i 



e-e' 



I d 3 fcp(fc)r; m (fc)^ m ,(fc)e i (k)er,(k). 

(9) 

The expression follows directly from Eq.([5]). One may 
arrive at the same result by inverting Eq.([2|) to obtain 
the ci£ m , expanding 0(k, e,rj) as a Legendre series and 
using the addition theorem (0} to eliminate the Legendre 
polynomials Pi when integrating over the direction in the 
sky. 
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It is useful introduce the spherical components as in jA = / dkk 2p( k ) r w eS,(fc)8„ A (fc) + u*G°,(k)G$(kj\ , 
lAckerman et all (|2007T > of the direction vector u J w L 1 v ' eK ' 1 v ' 1 WJ 

(14) 

and 



V2 



, n = n z , 



(10) 



since then one may write 



k-h = 1\ / — 
V 3 

We arrive at 



n+Y+^k) + n-Yf^ft) + nol 1 ? . (11) 



2i< 



itf = &^p(fc)| w |^ef(fc)e^(fc). (is) 



The first term in Eq.(JT2J) is the isotropic contribution. 
The second is the cross term, for which the geometric 
coefficients are given by 



Here the source integrals arc 



^ 12 ^ Wm' = n + C+ m , + n-Q Vm , + rioCww ( 16 ) 



h = / dkk 2 P(k) [@° e (k)] 2 , 



o 



(13) where 
I 



$£',£-1^ 
$t',t-l\ 



(1 + to-1)(1 + to) 
2(2*-l)(2*+l) £V 

f (l-m- l)(l-m) 
2(2*- 1)(2£ + 1) £V 

I (e-m)(e + m) 
'(2£-l)(2*+l) + ^ +li 





— TO H 


-1)(*- 


m + 2) 




2(2*- 


H)(2* 


+ 3) 


(J 


+ TO-f 


-!)(' + 


TO + 2) 




2(2*- 


H)(2* 


+ 3) 



(*- TO+ 1)(* + TO + 1) 

(2*+l)(2* + 3) 



(17) 
(18) 
(19) 



One can check that Q m -e'm' = Ci'm';£m- The last term 
in Eq. (fT2f is the autocorrelation of the anisotropic piece. 
The geomet ric coefficients 6w/"m' have been previously 
presented in lAckerman et all ( 20071 ). We have them with 
an extra minus sign for the off-diagonal component^. 
The factor also results in odd-parity correlations 

being imaginary. 

Although our cosmology features anisotropics, we as- 
sume that the underlying model is Gaussian. For a ped- 
agogic di scussion of these statistica l properties and their 
tests, see lAbramo k, Pereiral (|2010f) . 

3. ANISOTROPICALLY STRESSED DARK 
ENERGY 

The main reason for disregarding the anisotropic stress 
in the dark energy fluid might be that a minimally cou- 
pled scalar field, conventional parametrisation of the 
inflationary energy source, cannot generate anisotropic 
stresses. However, since there is no fundamental theo- 
retical model to describe dark energy, one might miss 
the behind physics of acceleration by sticking to the as- 
sumption of zero anisotropic stress. Such stresses are 
quite by viscous fluids, any high er spin fields and non - 
minimally coup le d scal ar fields toolApplebv et al.l (120101) 



iBattve fc Mbssl (120091): iRodriguesl (120081 1 : iCampanelli 
(I2009T) : iDimastrogiovanni et al l (12008): C oorav et al. 
(120081) : lAkarsu fc Kilind (120101) : ICooke fc Lvnden-Bel] 



1 The reason for this discrepancy was a forgotten i* factor in 
the Ackerman paper. 



(|2009|) . To study such a generic property with a many 
possible theoretical realizations, it is useful to employ a 
parametrisation of its physical consequences. 

An efficient way to describe possible deviations from 
perfect-fluid cosmology is to introduce the post-general 
relativity cosmologic al parameter w along the lines of 
I Caldwell et al.l (|2007t) . which is defined as the difference 
of gravitational potentials in the Newtonian gauge, 



t/j= (1 + w) 
where the line element reads 



ds 2 



a 2 {n) [-(1 + 2ij))dif + (1 - 2cj>)dx l dxi 



(20) 



(21) 



This parameter then appears as a cosmological general- 
ization of the post-Newtonian 7, for wh ich one has tight 
constraints from the Solar system scales IWiEl (|2001l ). An 
economic assumption is then that such a generalized pa- 
rameter depends, at all relevant scales, on the ratio of 
matter and dark energy densities, 



W = Wq- 



Pde 

PM 



(22) 



Then one has reasonable constraints on the wq from var- 
ious scales, and ~cuq of order one would imply a variety 
of phenomenology at different scales, ranging from Solar 
system physics to cosmology, just at the verge of detec- 
tion. Here we study on cosmo logical effects o f dark en- 
ergy and adopt the recipe 3 of iCaldwell et al.l (|2007l ) as 
the basis to parametrize the shear stress of dark energy. 

In particular, we will explore the case that the 
anisotropic stress has a preferred direction as in 
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iKoivisto fe Motal (12008a!) . For each Fourier mode of cos- 
mological perturbations, we write 



vj = i(k ■ h)zuo 



PDE 

pM ' 



(23) 



where n is the direction of the anisotropy and wq is 
real. Consider then the transfer function Eq.©. The 
0°(fc) would now be as usual. Thus it includes contribu- 
tions from both early and late universe. At the largest 
scales the Sachs- Wolfe effects are known to dominate the 
anisotropy sources. The second part would be given by 
the rotationally non-invariant part of the ISW contribu- 
tion, which in our present prescription is the following: 



Qf(k) = -iw Q J 



e -K(r,)_^_ ( PDE 



drj V Pm 



<t>k(ri)j 3i[kr(r))]dri, 
(24) 

where 4>k(v) is given by the standard computation. Wc 
just add the extra contributions due to Eq. (f24"|) to the 
sources from which to compute the correlators as de- 
scribed in the previous section. Therefore it becomes 
straightforward to determine the features in the CMB 
sky in this prescription. 

This parametrization describes a gradient-type mod- 
ification of th e effective CMB sources. We note that 
iTangenl (|2009f ) has determined the implication s of a 
super-horizon perturbation, and lErickcek et al.1 ()2008l ) 
considered such a spatial variation of the curvaton field 
at inflation. In our model, the anisotropy is formed dy- 
namically and becomes important with the dominance of 
dark energy. Explicitly, wc have a spontaneous modifi- 
cation of the effective CMB sources through the gradient 
operator as follows: 



V>(x) = [1 + (n • V)] 0(x) 



(25) 



when \n\ = zuo{pde/pm)- This amounts to shifting the 
Fourier modes of the perturbations exactly as prescribed 
in flUJ) and l[25jl. 



1 + i(k ■ n)wo 



PDE 
PD 



(26) 



Since the anisotropic part develops as a result of the evo- 
lution of the universe, it is a property of the transfer 
functions and not of the primordial spectrum of of per- 
turbations. One does not expect odd A£ couplings from 
primordial origin, since they violate parity. The presence 
of the imaginary unit is necessary for reality of the per- 
turbations, which can be checked as follows. Since the 
physical perturbation is a convolution of the primordial 
physical perturbation and the Fourier transformation of 
the transfer function, one notes that the latter should 
also be real. We get that 



^(x) 



Primordial 



x) 



I f°° k 
— 7? / dk— sin (kx) 
27r 3 7 x V ; 



1 + (x ■ n) [ cot (kx) 

kx 



(27) 



where <p(k) is the (real, isotropic) transfer function which 
depends only on the magnitude of the wavevector, and 
the right hand side is the anisotropic transfer function in 
the configuration space that retains it's reality. 



To recap our approach, wc introduced a mismatch of 
the two gravitational potentials in the Newtonian gauge. 
This mismatch, quantified by tu, is given by a gradi- 
ent along a preferred axis n. In the following, we will 
also allow scale dependence of this effect by introducing 
the spectral index nd - This parametrisation can then be 
used to constrain the presence of the such gradients in the 
late universe, since they would be seen in the CMB (prac- 
tically only) through their impact on the time-evolution 
of the gravitational potentials because of the ISW effect. 
Physically, these gradients could be caused by a large- 
scale inhomogeneity entering our horizon, spontaneous 
formation due to e.g. coherent magnetic fields or simply 
the possible imperfect nature of the dark energy field. 

To clarify the difference of our approach to all previous 
literature, let us mention that odd modulations may be 
considered to occur at three distinct levels. The temper- 
ature fi eld itself can be mo dulated, for a recent exam- 
ple see lAluri fe Jain! (|2011l) . This would effectively de- 
scribe some systematics in the data. Strangely enough, 
the primordial spect rum itself could contain parity vi- 
olating contribution IKoivisto fe Motal (|2011| ). That is 
consistent only in the context of noncommutative quan- 
tum field theory, and thus provides a unique signal for 
such high energy modific ations of the standard model 
iGroeneboom et all (1201 Of ). Finally, the cosmological 
structures may evolve statistically oddly, which is the 
case we focus upon here. 

4. MODEL FITTING 

In this section we describe in detail the confrontation 
of the model with the WMAP data. We will now discuss 
the method used to obtain the set of parameters which 
gives the best fit between our model and the observa- 
tions. Our basis is the evaluation of the Likelihood func- 
tion in a 4-dimensional parameter space. After explain- 
ing the general likelihood procedure, we explain in more 
detail the different steps taken to calculate and maximize 
the likelihood. 

4.1. Data model and notation 

Given a set of data {<ij our goal is to find the set 
of parameters which maximizes the posterior. For ease 
of notation let a = (9,<f>,mo,n<ie) — (n, vuq, n^c) denote 
the set of parameters to be determined. By Bayes' the- 
orem we know that the posterior distribution P(a|d) cx 
P(d|a)P(a) = £(a)P(a) where C{a) is the likelihood 
and P{a) is a prior. We take a conservative approach and 
assume that we know nothing prior about the anisotropic 
parameters, and thus P(a|d) = C(a) up to a normaliza- 
tion constant. If we manage to compute the likelihood 
function in all of the parameter space then we automat- 
ically have the posterior distribution and our job is es- 
sentially done. 

Although our model is anisotropic wc still assume that 
the underlying distribution is Gaussian, and thus the 
likelihood is 

D -id t C- 1 (a)d 



£(a) oc 



v /dctC(a) 



(28) 



where the data vector d consists of the ai m of the ob- 
served masked map. The correlation matrix C = S+N is 
the sum of the CMB signal covariancc matrix S and the 
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Fig. 1. — (Unnormalizcd) I ntegrals (c q, 1311 - 1331) computed from 
a modified version of CAMB ILewis et al, I 12000 ). Notice how the 
anisotropic integrals decay towards zero after only a few multipolcs. 
These integrals arc all used in the construction of the signal covari- 
ance matrix. An anisotropic scalar spectral index of n^e = 1-0 is 
used in this plot, and we have not normalized them. (Color version 
of this figure is available online) 

noise covariance matrix N. Our analysis is performed in 
harmonic space where the signal covariance Si m - t i> m > = 
{ae m a}, m ,} is computed from Eq. (fl"2|) and thus contains 
non-diagonal anisotropic contributions from dark energy. 
This matrix gives the dependence on the cosmological 
parameters. 

Finally the observed data vector d may be written in 
harmonics space as 

dim = h e w e s£ m + n fm (29) 

where hi is the instrumental beam, wi is the pixel win- 
dow function and n^ m is the (Gaussian) noise term. Since 
there is no correlation between the signal and the noise 
we have 

(d em d* e , m ,) = (§e m S* e , m ,) + (n im n* e , m ,) (30) 

where S£ rn = hi~wisi m is the observed signal. The goal 
is now to maximize this likelihood with respect to the 
model parameters, but we will first explain in some detail 
how we calculate the covariance matrices involved in the 
likelihood calculation. 

4.2. Signal covariance 

From equations (|12II19[) we see that the covariance ma- 
trix, in addition to the diagonal isotropic contribution Ii 
contains the cross term coefficients Qm-j'm' which couple 
I to £' = [£ ± 1] and m to to' = \m,m± I]. The last 
term ^em-j'm' is the Ackerman (jAckerman et al.l l2007| ) 
term for which we have couplings when £' = [£, £ ± 2] 
with ml = [to, to ± 1, to ± 2]. All other terms are zero. 

We will now express the integrals in equations (|17H19p 
in terms of power spectra Cg, C^,, Cf e ^ as 

Ci = -lu (31) 

C«' = — ifi' (32) 
Cff = (33) 



From equations (|17H19j) we see that Cf e , and Cf e ^ only 
give contributions for £' = I ± 1 and I' = £ ± 2. In 
figure [1] we have calculated (using a modified version of 
CAMB) and plotted these integrals (because of symme- 
try the £' = £ + 1 and £' — £ + 2 terms are equal to 
the £' = £ — 1 and t' —1—2 terms plotted) and com- 
pared to the isotropic power spectrum. From this fig- 
ure we clearly see that the anisotropic contribution from 
the dark energy component becomes negligible except 
at the largest scales where the anisotropic contribution 
even exceeds the isotropic. The most prominent alterca- 
tion comes from the anisotropic integral Cf^_ 2 contribut- 
ing to the off-diagonal elements of the signal covariance 
matrix. Due to the short range of the anisotropic in- 
tegrals we have altered the pivot scale in CAMB from 
k = 0.05 Mpc -1 to k = 2 x lO^Mpc -1 . In this way, 
we ensure that the spectral index enters correctly to tilt 
these integrals. 

4.3. Transformation of variables 

We have noticed from simulations that there was a 
significant degeneration between the anisotropic spectral 
index ride and the amplitude wq due to the fact that 
they both regulate the magnitude of the non-diagonal 
signal matrix elements. In order to ease the estimation 
procedure, we choose to estimate for a new variable Wq 
instead of zuq defined by w§ = zuq / \fa where we define 
a as a = A(I^ A 2 (1.0))/A(I^ A 2 (n dc )) where A(I*(n de )) 
is the area under the anisotropic integral Ijfe, (ride) where 
X = {A, AA} and ride is the dark energy spectral index. 
The parameters n^c and tuq are not degenerate and can 
thus more easily be estimated for. In the end, we convert 
to the physical parameter wq and all results are quoted 
in terms of this parameter. 

4.4. Modification of spectrum 

The angular power spectrum Ci will receive a contribu- 
tion from the anisotropy which could have an observable 
impact on the largest scales of the universe. This can be 
seen from equ ations (I17H19I) an d the f orm of the & m ;£'m' 
elements (see [Ackerma n et al.l (|2007l ) for details) . One 
must therefore be careful when performing the full anal- 
ysis, making sure that any choice for xuq does not signifi- 
cantly affect the power spectrum away from the WMAP 
best fit spectrum, but only the anisotropic contribution 
to the correlations between ai m s. 

To quantize these statements we calculate the net extra 
power from the anisotropic contribution. It is given by 
the diagonal part of the anisotropy which can be written 

2jAA w 2 1 

AC'l = " & m ;lm (34) 

m— — I 

Using the explicit form of t;i m -,im to perform the sum- 
mation we find the modified power spectrum (see the 
Appendix for details): 

Cr d = Ci + J-^4 A (n de ) = C t + |c, AA (35) 

The extra contribution to the power-spectrum from the 
anisotropic sources depends on the amplitude parame- 
ter (degree of isotropy breaking) and the dark energy 
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spectral index (parameterisation of the fluid scale- depen- 
dence) which we have included explicitly as an argument 

in T AA 
in i u . 

We see from equation 1551 that there is a limit to what 
values the amplitude zuq may take in order to obtain a 
power spectrum which is consistent with the WMAP7 
best fit. This is however only true when we assume the 
other cosmological parameters to have the WMAP best 
fit values. Clearly the new parameters which we have in- 
troduced allows for the other parameters to vary and one 
should rc-cstimate the other cosmological parameters to- 
gether with wq and ride- In order to test the anisotropic 
model without running a full cosmological parameter es- 
timation one may renormalize the covariance matrix for 
a given set of parameters (wq, n^e) in such a way that the 
power spectrum is kept constant at the best fit WMAP 
model. In this case the signal covariance becomes 



-iNorm 



s 



(36) 



J 



N, 



<imi;fcm2 



(-W (2<i+ 4r 2+i) x>-^ 

^3=0 



where St m lm — C e + ^m-jm^o^f ■ With this normal- 
ization, the diagonal part of our signal covariance matrix 
will match the WMAP power spectrum regardless of am- 
plitude and spectral index for the dark energy, while the 
off-diagonal components describe relative anisotropy. 



4.5. Noise covariance 
The noise in pixel space is assumed to be uncorrclatcd 



between pixels, i.e. N, 



{mrij 



of Sin where i and 



j are pixel indices, and is the noise root-mean-square 
deviation. The noise covariance matrix in pixel space is 
therefore diagonal. When going to spherical harmonic 
space, the harmonic coefficients of the noise are corre- 
lated and Ni m; e' m i = (ni m ng, m ,} is therefore a dense 
matrix. 

Expanding the noise harmonic coefficients in terms of 
pixel space quantities we eventually find that the expres- 
sion for the noise matrix in harmonic space becomes 



1 5(m 3 -m 1 +m 2 ) 

\ / \ 1Tl 3 —mi 7712 / 

(37) 



r 



where the Dirac delta-function arises due to the require- 
ment iiij — m 2 — rn 3 ~ for the Wigner 3j symbols. It is 
also required that the triangle condition \£i — £% \ < £3 < 
£1 + £ 2 is fulfilled. Notice also that due to this relation 
the sum over £ 3 goes up to 2£ max . 

The a^ 3 m 3 coefficients originate from a spherical trans- 
form of the variance of the noise map, of. Eq. (|5T|) is 
then implemented into our code. It needs only be com- 
puted once for each run of the code and is added to the 
signal covariance matrix in the step before the skycut is 
applied. 

4.6. Correlations introduced by the mask 
If we let Cim-e'm' denote the covariance matrix with- 
out mask, and Ci m -^' m ' denote the corresponding matrix 
including correlations from the sky cut, then the relation 
between them in harmonic space is found to be 

(38) 

LM L'M' 

which can be written compactly in matrix form as C = 
WCW'. The operation in Eq. (f3"5)) can be shown to be 
additive so the covariance matrix is the sum of the sig- 
nal plus noise correlation matrices. The multipole range 
here is L,L' £ [2,£ max ], and the sums over M,M' here 
run over positive values. The hermitean coupling matrix 
Wim-i'm' defined by 

wWm'= [ dnM(ii)Y? m (a)Yt, m ,(si) (39) 

is a function of the pixel space mask M(f2) (where f2 = 
{9, <j)) is the angular position on the sky) and so depends 
on the resolution iV s id c . It quantifies the new couplings 
between modes that arise due to the fact that we are now 
not analyzing a full sky. 



Starting with the WMAP KQ85 mask at A^ sidc = 512 
we degrade our mask so that the operation of applying 
the mask in pixel space can be traced exactly by applying 
the kernel matrix W in harmonic space. This is done by 
first smoothing with a Gaussian beam of fwhm= 744 ar- 
cmin, and then setting M(p) = (where p is a HEALPix 
pixel index) where M{p) < 0.80. We then band-limit 
the m ask so that it contains multipoles in the desired 
range (jArmendariz-Picon fc Pckowskv 2008h. These op- 
erations ensure that our mask does not contain small- 
scale structures. In the process the mask is expanded so 
that it now covers about 25% of the sky. 

It now remains to give an expression for the coupling 
kernel. Eq. (|39|) can be transformed by decomposing the 
mask into spherical harmonics and then performing the 
resulting integral over all angles to obtain again Wigner 
3j symbols. This is exactly the same analytical proce- 
dure which led to Eq. (|37[) with some minor modifica- 
tions. The result in this case becomes 

/ (2li + l)(2l a + l) 

Wt imi ;l 2 m 2 -K- 1 ) Y ^ 

2 ^ i / h h h \ ( h h ^ 

}2 a t3m3 y/2e 3 + l 

i 3 =o \ J \ m 3 mi -m 2 

5(m 3 +m 1 -m 2 ) (40) 

almost identical with Eq. (|37|) . As in the previous case 
the internal sum covers multipoles up to 2^ max , and the 
main difference is a change in sign of mi,77i2- We find 
that the mask coupling matrix is relatively well condi- 
tioned and further manipulations are not necessary. This 
has been tested by simply inverting the kernel matrix W 
constructed from the resulting mask which is used in the 
analysis of the WMAP data. 
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4.7. Likelihood maximization scheme 

In order to maximize the likelihood we use a non-linear 
Newton-Rapson search algorithm^ for the direction and 
amplitude. Finding the maximum of the likelihood is 
equivalent to finding the minimum of the quantity 

-21og£ = d t CT 1 d + TrlogC. (41) 

which has a global minimum when the likelihood func- 
tion has a global maximum. The algorithm minimizes 
a general unconstrained function by evaluating the first 

I 

0(-2iog£) = d jdc- 1 d + Tr /aiogC 

da da \ da 



and second derivatives. Due to the symmetry of the sig- 
nal covariance matrix, zuq is constrained to be larger or 
equal to zero. A negative amplitude can always be re- 
placed with a positive one and a shift in the angles (6, <fi): 
S(— wo, n^c, 0, (f>) = S(tuq, ride, 7r — n+cj)). This is easily 
seen from the definition of the signal covariance matrix. 

The gradient is computed analytically, the derivative 
of Eq. (|4ip with respect to any of the parameters in the 
set a is 



_ dtc -.§§ C-d + *(c->§§) (42) 



To find the derivative of C _1 we have differentiated the 
identity matrix I = CC 1 and solved for the derivative 
of the inverse covariance matrix in terms of the deriva- 
tive of the matrix itself. The analytic derivative of C 
is computed from Eq. [12J The second derivatives are 
estimated from the gradient using a secant method. De- 
pending on our initial guess for the parameters and our 
convergence criteria the minimizer in general may or may 
not converge to a global minimum. In our case false con- 
vergence is rarely a problem since the likelihood surface 
is well behaved and the local minimum has a small am- 
plitude compared to the global. 

For the spectral index ride we run a grid calculation. In 
each grid point, we apply the above maximization pro- 
cedure and find the value of the maximum likelihood for 
the given value of ride- In the end we search the grid to 
find the full global maximum in the 4-paramctcr space. 

5. APPLICATION TO WMAP-DATA 

Let us now discuss the results obtained with and with- 
out normalization of our signal covariance matrix (see 
section 14. 4p . We analyze the V-band (61 GHz) data 
map which is believed to be one of the cleanest bands 
in terms of foreground residuals, and recommended for 
cosmological analysis by the WMAP team. To this map 
we apply the modified WMAP KQ85 galactic skycut, re- 
moving 25% of the sky. Since we are only analyzing the 
largest scales no special care is taken with regards to 
point-source masking. We take into account the noise 
RMS pattern and the corresponding beam properties for 
the V-band. We analyze these maps out to a maximum 
multipole moment £ max = 20. The kernel matrices W 
which emulate the effect of a skycut in harmonic space 
includes multipoles up to £ max = 40. Now that we have 
our data map we are ready to start the analysis. 

5.1. Unnormalized covariance matrix 

Performing a grid-calculation with a spectral index 
range of — 5 < ride < 5 with a stcpsize of Aride = 0.1, 
where for each value of rid c we do a 3-dimcnsional search 
for the peak of the likelihood function using our likeli- 
hood maximization scheme described in section 14.71 we 
find that the likelihood for negative spectral index val- 
ues are very insignificant. As we approach ride = the 

2 using dmng.f from www.netlib.org 




2.5 3.0 3.5 

Spectral index n de 
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Fig. 2. — 2-dimcnsional plot of the raw likelihood (Posterior dis- 
tribution) as a function of dark energy spectral index n^ c and the 
transformed amplitude roJJ , using WMAP 7-yr data. The spherical 
angles have been fixed to their WMAP best-fit values at ^ max — 20 
to enable a projection. The posterior mode value is located at 
n dc = 3.1. The (transformed) amplitude found in this plot is 
higher than the value for the true anisotropic amplitude noted in 
table [T] due to a bias introduced by application of the mask and a 
linear transformation. 



likelihood starts peaking slowly until we find a peak at 
ride = 3.1. The best-fit direction remains practically 
constant as me move through the grid (the change is 
completely negligible compared to the uncertainty) indi- 
cating that correlation between the two sets (6*, <fi) and 
(ndc, wo) is weak. 

We find Fisher matrix error bars calculating the Fisher 
matrix using 



1. 



a/3 



-Tr 



^Lc- 1 — C- 1 
d\ a d\p 



(43) 



where the derivatives of the covariance matrix are found 
analytically for the direction and amplitude and numer- 
ically for the spectral index. Since the amplitude and 
spectral index are weakly correlated, the off-diagonal el- 
ements are taken into account in the matrix. The results 
are shown in table [T] As expected, in order for the model 
to be consistent with the power spectrum, we find vjq 
consistent with zero within the la error. In figure [2] we 
show the likelihood surface close to the peak. Note that 



Fig. 3. — Distribution of (0,<p) values on the sphere from Gaus- 
sian (too = input) simulations. Note how the positions found 
in the simulations are randomly distributed on the sphere and not 
aligned along some particular axis, in clear agreement with a ran- 
dom Gaussian distribution. 



TABLE 1 

Results from WMAP 7yr-data. Lower line shows the 

RESULTS WITH NORMALIZED COVARIANCE MATRIX. 



"■dc 


roo 


h(l,b) 


3.1 ± 1.5 


0.51 ±0.94 


(168°, -31°) 


1.2 ±0.7 


7.12 ±3.82 


(179°, -27°) 



the amplitude in this plot is the transformed amplitude 
(see section l4~3j) and that the amplitude at the maxi- 
mum of the likelihood is biased with respect to the best fit 
amplitude. This bias results from the complicated form 
of the likelihood introduced by the mask. The bias is 
corrected for in the following manner: given the param- 
eters found from the peak of the likelihood, we generate 
100 anisotropic realisations. For each realisation we esti- 
mate the anisotropic parameters and in the end compute 
the average bias in Wq. Next we subtract the bias from 
the input-value and repeat the procedure until our av- 
erage computed amplitude matches the value found in 
the WMAP data. When the bias has been subtracted 
we are left with "the true" estimate of the transformed 
amplitude Wq. The fiducial amplitude is then obtained 
using tuq = y/am^j. In the unnormalized case we find a 
final amplitude value zuq = 0.51. 

The best fit direction is somewhat close to the galactic 
center. In order to check that this is not caused by the 
shape of the mask, we estimated the direction on 1000 
simulated isotropic maps and found that there is no bias 
towards the galactic center. The estimated directions are 
shown in figure [3J 

The error in the amplitude has also been estimated 
using 1000 Monte Carlo simulations. We find that the 
error from simulations agrees well with the error found 
using the Fisher matrix. In figure 2] we show the best fit 
direction with error bars. 

An isotropic universe is clearly preferred by the data 
using this model. 

5.2. Normalized covariance matrix 

The model with an unnormalized matrix (see Eq |12[) 
giving the modified power spectrum in Eq. (|35|) is clearly 
not preferred by the WMAP 7-yr data. However, if we 
allow other cosmological parameters to vary we may be 
able to find a better fit as explained above. We have 
therefore repeated the procedure using the normaliza- 
tion in Eq. (|36|) which means fixing the diagonal part 



Fig. 4. — Map indicating the Icr uncertainty in the preferred 
direction of the axis. The background is the V-band (61 GHz) 
WMAP 7-yr data map with the KQ85 mask. The two axes plotted 
are the directions given in table [l] for the unnormalized and the 
normalized cases. 

of the covariance matrix to the best fit Cg regardless of 
amplitude and spectral index. 

The lower line in table [1] shows the results for the 
anisotropic cosmological parameters from the explo- 
ration of the likelihood space with a normalized covari- 
ance matrix. Again, a grid was set up for the spectral 
index in the interval — 5 < n^c < 5 with the same step- 
size and for each value of n^c we estimated the best-fit 
values of (9,<p,vjo) and the corresponding value of the 
likelihood. The preferred amplitude in the normalized 
case is so large that a huge change of cosmological pa- 
rameter values would be necessary to obtain the best fit 
WMAP spectrum. No physical solution is found in the 
normalized case. 

6. CONCLUSIONS 

In this work we tested anisotropic dark energy models 
with the 7- year WMAP temperature observations data. 
If dark energy is not a perfect fluid but for instance a vec- 
tor field, the CMB sky will be distorted anisotropically 
on its way to us by the ISW effect. The signal covariance 
matrix then becomes nondiagonal for small multipolcs, 
but at t > 20 the anisotropy is negligible. This can be 
used to constrain violations of rotational invariance in 
the late universe, and to obtain hints on possible imper- 
fect nature of dark energy and the large-angle anomalous 
features in the CMB. 

To model this phenomenon, we introduced a mismatch 
of the two gravitational potentials in the Newtonian 
gauge. The mismatch, quantified by w, is proportional 
to a gradient along the preferred axis n. We also al- 
lowed this effect to depend on the scale by introducing 
the spectral index n^c- Physically, such gradient could be 
caused by a large-scale inhomogeneity entering our hori- 
zon, spontaneous formation due to e.g. coherent mag- 
netic fields or simply the possible imperfect nature of 
the dark energy field. Many possible realisations of the 
latter possibility were discussed in the introduction and 
in the section [3] The dominant effect on the CMB is 
then a quadropole modulation, which has the same geo- 
metrical correlation structure but different time and scale 
dependence than in models considered previously. Now 
a dipolc modulation, though subdominant, is predicted 
too. 

We calculate the mode couplings introduced to the 
spherical harmonic coefficients of the CMB by the 



anisotropic model and obtain the full likelihood for the 
lowest multipoles where the dominant contribution to the 
model is expected to be found. Maximizing the likeli- 
hood taking into account the instrumental parameters 
of the WMAP experiment, we are able to find optimal 
estimates of the anisotropic parameters. Analysis of the 
masked WMAP V-band, fixing other cosmological pa- 
rameters, gave a best fit amplitude vuo = 0.51 ± 0.94 and 
"■do = 3.1 consistent with an isotropic universe. 

In comparison, test of the isotropic version of this 
parametrisation show that the data is then compatible 
with a vanishing de viation, and a l lows a nonzero vo of 
the order of O(0. 1) (|Daniel et all ()2010D ). At the level 
of Solar system, no hints of deviations are observed, 
and the post-New tonian corre ction is constrained to be 
at most O(10 -5 ) (jWilll ()2001h ). However, the numbers 
themselves cannot be directly compared, since our best- 
fit model features also a strong scale-dependence of the 
deviation. The preferred steep blue spectral index may 
be due to the fact that the statistically most significant 
contribution must come from smaller scales, the largest 
scales of the CMB being severely cosmic variance limited. 

Another shortcoming of our parameterization is its in- 
ability to incorporate the lack of large-angle power in the 
observed sky, one of the most striking anomaly present 
in the data. Indeed, the main lesson to be derived from 
our study is that cosmological constraints on realistic 
imperfect post-Newtonian deviations are on the order 
0(1O -4 ) (reported in table [1]), only an order of magni- 
tude below those obtained from the Solar system scales. 
This motivates to further investigate the possible ori- 
gin and constraints of imperfect source terms in cosmol- 
ogy. In particular, a fully consistent post-Friedmannian 
param etrisation alo ng th e lines of iFerreira fc SkordTsI 
(|2010( ); lBaker et al.1 (|2011l ). tailored to the study of direc- 
tional dependence of deviations from the standard pre- 
dictions of linearised cosmology, remains to be developed. 
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APPENDIX 

ANISOTROPIC CONTRIBUTION TO THE POWER SPECTRUM 

The diagonal part of the covariance matrix is a sum of the power spectrum due to the isotropy and a term determined 
by the anisotropic parameters mo, nde and : 

2 

The dependence on the spectral index comes from the integr al over the anisotropic transfer functions in i^ A . The 
diagonal part of the geometric factor is ([Ackerman et a l. 2007) 

-I + £(£ + 1) + m 2 2 2£(£ + 1) - 2m 2 - 1 
W™- Ai+n_ (2£ _ 1)(2£ + 3) + n o {2 £-l)(2£ + 3) (A2j 



where the spherical components n + ,n_,no containing the angular dependence have been defined in Equation (|10p . 
Using the well-known result 

" 2 £(£+l)(2£+l) 
2^m = (A3) 

m=l 

we find that the average of the geometric factor is free from angular dependence and simplifies nicely to 

2 / «m;«m — ~ (A4) 

m=-i 

With this result one finds that the theoretical prediction for the modified power spectrum due to the anisotropic 
component becomes 

CT ° d = 2fTi E< *»°*"> = ^ + t qAA (A5) 

m 

which is the modified power spectrum C" lod quoted in Eq. ([35]) . 
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